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Abstract Bond-percolation graphs are random subgraphs of the d-dimensional in- 
teger lattice generated by a standard bond-percolation process. The associated graph 
Laplacians, subject to Dirichlet or Neumann conditions at cluster boundaries, repre- 
sent bounded, self-adjoint, ergodic random operators with off-diagonal disorder. They 
possess almost surely the non-random spectrum [0, Ad] and a self-averaging integrated 
density of states. The integrated density of states is shown to exhibit Lifshits tails at 
both spectral edges in the non-percolating phase. While the characteristic exponent of 
the Lifshits tail for the Dirichlet (Neumann) Laplacian at the lower (upper) spectral 
edge equals d/2, and thus depends on the spatial dimension, this is not the case at the 
upper (lower) spectral edge, where the exponent equals 1/2. 
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Introduction 

Spectral graph theory studies linear operators which are associated with graphs. The 
goal is to see how properties of the graph are reflected in properties of the operators 
and vice versa. This has attracted vivid interest in the last two decades [21,15,12, 
14]. The kind of graphs we shall be concerned with in this paper are bond-percolation 
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graphs [16], a special type of random subgraphs of the d-dimensional integer lat- 
tice. They are of popular use in Physics for modelling various types of random en- 
vironments [28,9]. On these graphs we consider Laplacians with different kinds of 
boundary conditions at cluster borders and study some of their spectral properties. 
Apart from the non-randomness of the spectrum, and of the spectral components in 
the Lebesgue decomposition, we show the existence and self-averaging of the inte- 
grated density of states. The main result establishes Lifshits tails of the integrated 
density of states at the lower and upper spectral edge in the non-percolating phase. 
Depending on the boundary condition and on the spectral edge, the Lifshits tail dis- 
criminates between stretched (i.e. linear) and condensed (i.e. cube- or ball-like) clus- 
ters which contribute the dominating eigenvalues. The crucial technical estimates in 
our proof are Cheeger [14] and Faber-Krahn [13] isoperimetric inequalities on graphs. 
Our analysis here is facilitated by the fact that in the non-percolating phase almost all 
graphs consist of infinitely many finite clusters. Yet, the non-percolating phase gives 
rise to interesting phenomena, because it supplies clusters of arbitrarily large size. In 
the percolating phase one has to cope with the infinite cluster, too, which requires a 
more intricate understanding. This case will be studied in [23]. Related previous work 
can be found in [1,5]. 

Spectral properties of Laplacians on bond-percolation (or related) graphs have 
been studied in the Physics literature, see the general accounts [28,9] or the recent 
examples [4,8,22] for applications to soft matter. The Lifshits tails, whose existence 
we prove here, were sought after in the numerical simulations [7] for the Neumann 
Laplacian on two- and three-dimensional bond-percolation graphs. The tails could 
not be observed there due to finite-size corrections and the considerable numerical 
effort needed to access such rare events. For the different case of Erdos-Renyi ran- 
dom graphs, however, the existence of Lifshits tails for the Neumann Laplacian was 
known on the basis of analytical, non-rigorous arguments [6,7], which inspired our 
proof here. Moments of the eigenvalue density for this model were rigorously anal- 
ysed in [17]. Other models in the physics literature deal with the adjacency operator 
on bond-percolation graphs [18,26]. Quite often, this goes under the name quantum 
percolation. Yet, from a rigorous mathematical point of view, all of the above models 
with off-diagonal disorder have remained widely unexplored, see however [1]. 

In contrast, Laplacians on s/fe-percolation graphs of the rf-dimensional integer lat- 
tice belong to the class of models with diagonal disorder. Therefore they are closer 
to the range of applicability of the highly developed theory of random Schrodinger 
operators [19, 10,24,29,20]. This is partly of help for analysing their spectral proper- 
ties with mathematical rigour. For finite-range hopping operators on site-percolation 
graphs, the non-randomness of the spectrum and existence of the integrated density 
of states was shown in [30,31]. Particular emphasis is laid on the behaviour of the 
spectrum related to finitely supported eigenfunctions, see also [11], where the issue 
was first taken up from a mathematical point of view. Furthermore, a Wegner estimate 
is established in [31] for an Anderson model on site-percolation graphs. Highly de- 
veloped large-deviation techniques for the parabolic Anderson model are used in [5] 
to prove Lifshits tails for the integrated density of states of the Laplacian A g (in the 
sense of our Definition in Eq. (1.9) below) on s/fe-percolation graphs. 
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This paper is organised as follows. In Section 1 we give the precise definitions of 
the objects we are dealing with and state our results. Theorem 1.14 in Subsection 1.3 
contains the central result on Lifshits tails of the integrated density of states. All proofs 
are deferred to Section 2. 



1 Definitions and Results 

1.1 Bond-percolation graphs 

A thorough and comprehensive account of (bond) percolation can be found in Grim- 
mett's textbook [16], which serves as a standard reference on the subject. For d G N , 
a natural number, we denote by h d the (simple hypercubic) lattice in d dimensions. 
Being a graph, the lattice L d = (Z d ,E d ) has a vertex set, which consists of the d- 
dimensional integer numbers Z d , and an edge set E d given by all unordered pairs 
[x, y] of nearest-neighbour vertices x, y G Z d , that is, those vertices which have Eu- 
clidean distance \x — y\ := (Y^t=i \ x v ~ vA 2 ) 1 ^ = 1- Here, elements of Z d are 
canonically represented as d-tuples x — (x\ , . . . , xj) with entries from Z. 

Given any subset of vertices / f C Z 11 and a subset of edges S C {[x, y] G 
E d : x,y e Y} between them, we call the graph := (Y, g) a subgraph of L d . The 
vertex degree 

%f{x) :=\{y&1 d : [x,y] £ g}\ (1.1) 

of x G Jj d counts the number of edges in & that share the vertex x as an endpoint. 
Here, \A\ denotes the cardinality of a subset A C Z d , and we use the convention 
|0| = for the empty set. A graph is called finite, if \Y\ < oo. 

A given graph consists of finitely or infinitely many clusters ^j, j = 
1, 2, . . . , J < oo, which are the maximally connected subgraphs of 'S . More precisely, 
tfj : = (fj, S'j) is a connected subgraph of Sf, if ^ C f,^ C g and if for every 
pair x, y £ Yj with x ^ y there exists K G N and Xk G Yj, k = 0,1,2, ... , K, such 
that Xo := x, xk '■= y and [xk-i,, Xk} E gj for all fc = 1, 2, . . . , K. A connected 
subgraph ^ of is maximal, and hence a cluster, if for every connected subgraph 
c g' := (Y', g') of & obeying Y' D Yj and g' D gj one has T = This defini- 
tion of clusters also includes isolated vertices, in which case ^ = ({x}, 0) for the 
corresponding x 6 Y with dy (x) = 0. Apparently, the decomposition of ^ into its 
clusters is unique - apart from enumeration. 

Next, we consider the probability space fi = {0, 1} E , which is endowed with 
the usual product sigma-algebra, generated by finite cylinder sets, and equipped with 
a product probability measure P. Elementary events in Q are sequences of the form 
ui = (w[ x ,j/])[x,i/]eE <i > an d we assume their entries to be independently and identically 
distributed according to a Bernoulli law 

F(w [XtV] = l)=p (1.2) 

with parameter p g]0, 1[, the bond probability. To a given ui G fl we associate the 
edge set 

g^ :— {[x,y] G E d : u>[ x , y ] = l} • (1.3) 
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Definition 1.1. The mapping Q 3 uj ^ Sf(^) := (Z d , S^) with values in the set of 
subgraphs ofh d is called bond-percolation graph in Z d . 

The most basic properties of bond-percolation graphs are recalled in 

Proposition 1.2. For d ^ 2 there exists p c e]0, 1[, depending on d, such that 

(i) for every p e]0,p c [, the non-percolating phase, one has 

pLert: ^ M consists of oo-many 1 =1 (14) 
[ clusters, which are all finite J 

(ii) for every p E]p c , 1[, the percolating phase, one has 

^ g q . consists of exactly one infinite 1 _^ ^ ^ 

cluster and oo-many finite clusters J 

Remarks 1.3. (i) The proposition collects results from Thms. 1.10, 1.11,4.2 

and 8.1 in [16], which were mainly obtained by Hammersley in the late fifties. The 
uniqueness of the infinite cluster, however, was only proven thirty years later by Aizen- 
man, Kesten and Newman, see [16]. 

(ii) In the one-dimensional situation, d = 1, one has p c = 1 and part (i) of the 
proposition remains true. 



1.2 Graph Laplacians 

The subsequent definition introduces Laplacian-type operators associated with an ar- 
bitrary subgraph of the integer lattice. The particularisation to operators on bond- 
percolation graphs follows at the end of this subsection. 

For a given subset A C Z d let I 2 (A) denote the Hilbert space of complex-valued, 
square-summable sequences that are indexed by A. 

Definition 1.4. Given any subgraph = (Y, <§) of h d with Y ^ 0, we introduce 
the following bounded and self-adjoint linear operators on £ 2 (Y). 

(i) The degree operator DifS) is defined as the multiplication operator with 
the vertex-degree function dcg : V — > N U {0}, x d^{x), that is, 

[D(&)<p}(x):=d w (x)<p(x) (1.6) 

for all ip e I 2 (Y ) and all x E Y. 

(ii) The adjacency operator A(&) is defined through its action 

for all ip e £ 2 (Y) and all x 6 Y. Here, we use the convention ^2 ye0 <p{y) = Ofor 
the empty sum. 
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(iii) The Neumann Laplacian is defined by 

Ax{9) := D(Sf) - A(Sf) . (1.8) 

(iv) The Pseudo-Dirichlet Laplacian is defined by 

Z\ B (Sf) := (2dl - D(Sf)) + 4 N (Sf) = 2dl - A(Sf) , (1.9) 

where t = ly stands for the identity operator on £ 2 (Y). 

(v) 77*e Dirichlet Laplacian is defined by 

A D (9) ■= 2(2dl - D(Sf)) + 4 N (SP) = 2dl + (2dl - D(Sf)) - A(Sf). (1.10) 

Remarks 1.5. (i) The asserted boundedness and self-adjointness of the Lapla- 

cians Ax(S^), X £ {N, D, D}, follow from the corresponding properties of D(@) and 
A(&). Indeed, since < dy{x) < 2d for all x £ V, it is clear that is self- 

adjoint and obeys < D(&) ^ 2dl in the sense of quadratic forms. The operator 
A(&) is symmetric, because 

<V,A(3»=2 1>*W<P(V) (1-11) 

[x,y]€g 

for all ip,ip £ £ 2 (Y), where (-0, <p) := J2 xe y ^*i x ) f( x ) denotes the standard 
Hilbert-space scalar product on £ 2 (Y). The factor 2 in (1.11) reflects that the sum is 
over unordered pairs. Moreover, applying the Cauchy-Schwarz inequality to (1.11), 
yields the upper bound 2d for the usual operator norm of A(&), and self-adjointness 
follows from symmetry and boundedness. 

(ii) The Neumann Laplacian Z\ N (S?) is called graph Laplacian or combinato- 
rial Laplacian in spectral graph theory, where it is commonly studied in various forms 
[21,15,12,14]. 

(iii) The quadratic form 

(<p,A N (9)<p)= J2 - f(y)\ 2 , <pe£ 2 (Y), (1.12) 



for the Neumann Laplacian reveals that 

Z\n(S?)^0. (1.13) 

Thus, a necessary and sufficient condition for ip £ £ 2 (Y) to belong to the zero- 
eigenspace of A N (%?) is that ip stays constant within each of the finite clusters of 
& (separately). Consequently, each finite cluster of & contributes exactly one zero 
eigenvalue to A^(^). In contrast, zero is not an eigenvalue of A^(^) and Ad(@). 

(iv) Let X £ {N, D, D} and let tfj := (Yj, j = 1, . . . J < oo denote the 
clusters a graph Sf := (Y, <§) is composed of. Then Ax(&) is block-diagonal with 
respect to the clusters, 

A x (9) = A x (tfj) (1.14) 



6 



Werner Kirsch , Peter Muller 



on £ 2 (y). We note that if ^ corresponds to an isolated vertex, then AyS^f) acts 
as multiplication by 7x on the one-dimensional subspace ^ 2 (^), where 7n := 0, 
75 := 2d, respectively 70 := 4d. 

(v) The Neumann and the Dirichlet Laplacian are related to each other. To see 
this we define a unitary involution U = U^ 1 = U* on £ 2 (Y) by setting 

(U(f)(x) := (-l)£2=i \ x »\tp(x) (1.15) 

for all a; G ^ and all ip e i 2 {"V). This involution commutes with and anti- 

commutes with A{&) so that D{&) = UD(^)U and 

A(&) = -UA(&)U. (1.16) 

Hence, we infer the relation 

Z\ D (Sf) =4dl -Z7Z\ N (S?)[7. (1.17) 
Combining (1.17) with (1.13), (1.9) and (1.10), we arrive at the chain of inequalities 

0^A N (&)^A b (&)^A D (&)^4dl. (1.18) 

(vi) Our terminology of the Laplacians is motivated by Simon [27]. Divide a 
graph & = (~V,g) into two subgraphs & k = (%, g k ), k = 1, 2, such that = 
and fi U ^2 = ^- Then one gets super-, respectively subadditive behaviour 

A N (9) ^ Z\n(^i)0^n(%), 

(1.19) 

z\ D (^)<z\ D (^)ez\ D (%) 

as a consequence of (1.12) and (1.17). Thus, introducing a separating boundary 
surface lowers Neumann eigenvalues and raises Dirichlet eigenvalues - in anal- 
ogy to the well-known behaviour of Laplacian eigenvalues of regions in continuous 
space, see e.g. Prop. 4 in Chap. XIII. 15 of [25]. In contrast, the eigenvalues of the 
Pseudo-Dirichlet Laplacian behave indifferently with respect to this proce- 

dure. Though, Z\g(Sf) is commonly termed a Dirichlet Laplacian in the literature. 

Next, we associate Laplacians to the bond-percolation graphs of Definition 1.1. 

Definition 1.6. The mapping Ax ■ ft 3 w <— > A^ := Ax(&^) with values in the 
bounded, self-adjoint operators on l 2 {T, d ) is called Neumann or Pseudo-Dirichlet or 
Dirichlet Laplacian on bond-percolation graphs in Z d , depending on whether X stands 
for N or D or D. 
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1.3 Results 

To begin with we summarise the most basic spectral properties of the Laplacian on 
bond-percolation graphs in 

Lemma 1.7. Fix X G {N, D, D} andp G]0, 1[. 7%en 

(i) the random operator Ax is ergodic with respect to Z d -translations. 

(ii) its spectrum is F-almost surely non random, more precisely, it is given by 
spec(Zix) = [0,4d] f-almost surely. 

(iii) the components in the Lebesgue decomposition of the spectrum are also 
¥ -almost surely non random. For every k G {pp, sc, ac} there exists a closed subset 

C R such that spec^Ax) = ^x^ f-almost surely. 

(iv) in the non-percolating phase, p G]0,p c [, the spectrum of Ax is F-almost 
surely only a dense pure-point spectrum with infinitely degenerate eigenvalues. The 
dense set of eigenvalues is also non random F-almost surely. 

Remarks 1.8. (i) The lemma is proven in Section 2. 

(ii) Part (ii) implies that the discrete spectrum of Ax is P-almost surely empty. 

(iii) As compared to the non-percolating phase considered in part (iv), there 
are additional spectral contributions from the percolating cluster if p G]p c , 1[. Among 
others, the percolating cluster contributes also infinitely degenerate, P-almost surely 
non-random eigenvalues corresponding to compactly supported eigenfunctions. This 
can be established with the same mirror techniques as it was done for related models 
on site-percolation graphs [ 1 1 , 30, 3 1 ] . Non-rigorous arguments [ 1 8, 26, 4] suggest the 
existence of continuous spectrum if p lies above the "quantum-percolation threshold" 

Pq > Pc- 

We proceed with the existence and self-averaging of the integrated density of states 
of Ax- To this end let S x £ £ 2 (Z d ) be the sequence which is concentrated at the point 
x G Z d , i.e. S x (x) := 1 and S x (y) :— for all y ^ x G T, d . Moreover, stands for the 
Heaviside unit-step function, which we choose to be right continuous, viz. 0(E) := 
for all real E < and 0(E) := 1 for all real E ^ 0. 

Definition 1.9. For every p G]0, 1 [ and every X G {N, D, D} we call the function 

N x : R 3 E hh. N X (E) := [ P(dw) {So, 0{E - A^ } )5o) (1.20) 

JQ 

with values in the interval [0, 1] the integrated density of states of Ax- 
Remarks 1.10. (i) Thanks to the ergodicity of Ax with respect to 7L d - 
translations, one can replace <5 n by 5 X with some arbitrary x G 7L d in Definition 1 .9 
without changing the result. 

(ii) The integrated density of states iVx is the right-continuous distribution 
function of a probability measure on R. The set of its growth points coincides with 
the P-almost-sure spectrum [0, Ad] of Ax- 



8 



Werner Kirsch , Peter Muller 



(iii) Forp < p c the growth points of Ax form a dense countable set, where Ax 
is discontinuous. These jumps in Ax are due to the infinitely degenerate eigenvalues 
of Ax, which arise solely from the finite clusters, cf. Lemma 1.7(iv). For p > p c there 
are also contributions to the jumps that arise from the percolating cluster. In addition, 
the set of growth points of Ax should not be restricted to discontinuities for p > p c , 
cf. Remark 1.8(iii). 

(iv) Eqs. (1.16) and (1.17) imply the symmetries 

N b (E) = l- lim N b (e), 

sHd - E (1.21) 
Nv W (E) = l-Jm_ E N ni > ) (s) 

of the integrated densities of states for all EeR. 

Definition 1 .9 of the integrated density of states coincides with the usual one in terms 
of a macroscopic limit. To make this statement precise, we have to introduce restric- 
tions of Ax to finite volume. 

Definition 1.11. Let & = (Y, S)be a subgraph ofl/ and consider a subset A C 1 d . 

(i) The graph <g A := (Y A , S A ) with f A := V n A and S A := {[a;, y] G £ : 
x, y G ~f A } is called the restriction of Sf to A. In particular, = (A, is the 
restriction to A of a realisation = (Z d , S^) of the bond-percolation graph. 

(ii) For X G {N, D, D} we define the restriction Ax yA of the Laplacian Ax to 
( 2 {A) as the random operator with realisations A^ A '■= Ax{y ( A ] )foralluj G Q. 

Lemma 1.12. Given p g]0, 1[ and X G {N, D, D}, there exists a set Q' C Q of full 
probability, P(j?') = 1, such that 



Nx(E) = lim 



|1 t^ce e2{A) 0(E-A^ A ) 



(1.22) 



holds for all lo G fi' and all E G M, except for the (at most countably many) discon- 
tinuity points ofNx- 

Remarks 1.13. (i) As to the limit A j Z d , we think of a sequence of cubes 

centred at the origin whose edge lengths tend to infinity. But there exist more general 
sequences of expanding regions in Z d for which the lemma remains true. 

(ii) The proof of Lemma 1.12 for X G {N, D} follows from the Ackoglu- 
Krengel superergodic theorem on account of (1.19), see Thm. VI. 1.7 in [10] and the 
discussion after Eq. (VI. 16) there. For X = D, the proof follows from Lemma 4.5 in 
[24], which establishes weak convergence of the associated density-of-states proba- 
bility measures, and Thm. 30.13 in [2]. 

(iii) The arguments in Sec. 6 of [31] show that the convergence in (1.22) holds 
whenever E is an algebraic number, that is the root of a polynomial with integer 
coefficients. Hence, the convergence (1.22) may even hold at discontinuity points of 
A^ X - In particular, forp < p c it holds for all E eR. 
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The central result of this paper is 

Theorem 1.14. Let d G N and assume p G]0,p c [. Then the integrated density of states 
Nx of the Laplacian Ax on bond-percolation graphs in 7L d exhibits Lif shits tails at 
both the lower spectral edge 

,. ln\\n[N x (E) - N X (0)}\ / -1/2 X = N, 

!™ — 1 r-F 1 = i a In f° r v ™ ( L23 ) 



E|o InS \-<i/2 ' X = D,D 

and at the upper spectral edge 

,. ln|ln[Ar-(4rf)-iVx(£)]| / -1/2 X = D, 

J?& M4^£) = l^/2 ^ X = N,D, (L24) 

w/zere N^{M) := \mi E ^ id N x {E). 

Remarks 1.15. (i) The Lifshits tails at the upper spectral edge are related to 

the ones at the lower spectral edge due to the symmetries (1.21). 

(ii) Remark 1.5(iii), the symmetries (1.21), Lemma 1.12 and Remark 1.13(iii) 
imply the values 

N b (0) = N u (0) - 0, iV N (4d) = N b (4d) = 1, 

1 - ^ (4d) = AVO) = Hm ^ = «(p) 

/tTZ d 

for the constants in Theorem 1.14. Here, n(p) is the mean number density of clusters, 
see e.g. Chap. 4 in [16]. Thanks to the right-continuity of the Heaviside function, the 
operator Q(—A^\) is nothing but the projector onto the null space of A^\. 

(iii) The Lifshits tails for N-^ at the lower spectral edge - and hence the one 
for JVd at the upper spectral edge - is determined by the linear clusters of bond- 
percolation graphs. This explains why the associated Lifshits exponent —1/2 is not 
affected by the spatial dimension d. Technically, this relies on a Cheeger inequality 
[14] for the second-lowest Neumann eigenvalue of a connected graph. 

(iv) If d ^ 2, then all other Lifshits tails of the theorem are determined by 
the most condensed clusters of bond-percolation graphs, like cubic clusters (see Re- 
mark 2.5 below for their definition), as they maximise the mean vertex degree among 
all clusters with a given number of vertices. In the proof of the theorem this will follow 
from a Faber-Krahn inequality [13] for the lowest (Pseudo-) Dirichlet eigenvalue of a 
connected graph. In contrast, for d = 1 there are no other clusters than linear ones, 
and the Lifshits exponent cannot discriminate between different boundary conditions. 

(v) For^/fe-percolation graphs, a stronger statement than (1.23) is known for 
the case X = D, see [5]. 
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2 Proofs 

In this section we shall prove Lemma 1.7 and Theorem 1.14. 

Proof (of Lemma 1.7). We follow the standard arguments laid down in [19,10,24]. 
The function fi 3 uj i— > A^\ which takes on values in the set of bounded self- 
adjoint operators on £ 2 (Z d ), is measurable, and the probability measure P is er- 
godic with respect to the group of translations (r z )^ eZ d on fi, which act as t z oj := 
i UJ [x+z.y+z])[x.y]ez d - Moreover, for every z G Z d let T z be the unitary translation op- 
erator on l 2 (Z d ), that is, T z (p(x) := <p(x - z) for all ip G f 2 (Z d ) and all x G Z d . 
The operator identity A%" u) = T^A^T; holds for all z G Z d and all u G fi 
and renders Ax an ergodic random operator [19], as claimed in part (i) of the lemma. 
Part (iii) is now a consequence of the general theory of ergodic random operators [19, 
10,24]. 

As to part (ii) it suffices to show the inclusion 

spec(Z\^ } ) D [0,4rf] for P-almost every w G fi (2.1) 

for all X G {N, D, D}, because the opposite inclusion is already supplied by (1.18). 
To verify (2.1), we define the event 

fi := |w G fi : for every I G N there exists a cube C Z d 

with l d points such that S? ( ?2, = ) . (2.2) 

A i A i > 

Here, we say that a subset of Z d is a cube with l d points (or, equivalently, with edges 
of length / — 1 G N), if this subset is some translate of the <i-fold Cartesian product 
{1, . . . , l} d . Colloquially speaking, the condition in (2.2) requires all bonds inside of 
to be present. 

Now, fix an arbitrary E G [0, 4d] = spec(Z\) in the spectrum of the ordinary lattice 
Laplacian A = zix(L d ). Then, there exists a Weyl sequence (V'_E,n)neN C £ 2 (Z d ) 
for A, that is, ||^B,n|| == {i>E, n , ipE.n) 1/2 = 1 for all n G N and 

lim 11(^-^1)^11=0. (2.3) 

We may also assume without loss of generality that the support supp ipE,n is compact 
for all n G N, since A is bounded. Furthermore, if (ipE.n)neN is suc h a Weyl sequence, 
then so is [T Zn ^pE,n)nen with arbitrary z n G Z d . Thus, given any u) G fi there exists 
a Weyl sequence (V'^'^neN f° r A with the property that, loosely speaking, all the 
supports are contained well inside the cubes of (2.2). More precisely, we mean that 
given every lo G fi and every n G N there must exist an integer > 3 and a cube 
A ( ^\ from (2.2) such that min{ \x - y\ : x G supp^^, y G Z d \ A^ } > 1. This 
yields 

\\(A^-Et)^J = 11(4-^1)^11 (2.4) 
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for all n E N and all u E fl. Hence, (ip^ n ) n efi is a l so a Weyl sequence for A^ and 
we have shown the inclusion in (2.1) for all w E ft. But P(J?) = 1, as we shall argue 
now. 

For every given integer I > 2 let (Ai^)^^ C Z d be a sequence of cubes in 
1 d with l d points such that Ai, fJil n Ai^ 2 = 0, whenever [i\ ^ /i2- Then, the 
events Qi,^ := {u> E J? : = L^, } are pairwise statistically independent, 

and P(Hi^) > does not depend on \i E N. So the Borel-Cantelli lemma implies 
P(i7;) = 1 for all integers I ^ 2, where fl\ := limsup^^ Qi^. The proof of part (ii) 
is completed by noting that fl D Hi-i^Hi. 

Finally, we turn to part (iv) and assume p < p c . We observe that, in P-almost 
every realisation of a bond-percolation graph, the translates of any given finite cluster 
occur infinitely often. This follows from a Borel-Cantelli argument like the one in 
the previous paragraph. Hence, the block-diagonal structure (1.14) of Ax implies that 
the set of eigenvalues of Ax is P-almost surely given by the union of the spectra of 
Axi^), where ^ runs through all possible finite clusters in L d . In particular, the set 
of eigenvalues is a non-random dense set and all eigenvalues are infinitely degenerate. 
□ 

The remaining part of this section concerns the proof of Theorem 1.14. It re- 
lies on deterministic upper and lower bounds for small eigenvalues of clusters. The 
lower bounds are discrete versions of well-known isoperimetric estimates for Lapla- 
cian eigenvalues on manifolds. 

Definition 2.1. For X E {N, D, D} and a connected subgraph := (Y, S) o/L d 
with \y\ > 2 vertices, let denote the lowest non-zero eigenvalue of Ax{^)- 

Proposition 2.2 (Cheeger inequality). Let ^ := (Y, <§) be a connected finite sub- 
graph ofh d with \Y\ ^2 vertices. Then its lowest non-zero Neumann eigenvalue 
obeys 

4%) > [^H!. (2 5) 

The quantity h ch (&) :— min^-- \d-Jtf\/\J4f\ is the Cheeger constant, where the min- 
imum is taken over all subgraphs Jtf of$ whose vertex set W obeys \W\ ^ \Y\/2. 
Here, \dJff \ := {[x, y] £ § : x E W, y E Y \ W} denotes the edge boundary of ,yf 
in<S. 

Remarks 2.3. (i) Proposition 2.2 just quotes a special case of a more general, 

well-known result in graph theory, see e.g. Thm. 3. 1(2) in [14]. 

(ii) The simple lower bound ft-ch(^) ^ l/(\Y\/2) on the Cheeger constant 
yields 

4 X) W^^- (2-6) 

This bound produces asymptotically the correct \V\ -dependence as \f\ — > oo, if Sf is 
a linear cluster Jz?„, i.e. a connected subgraph of h d having 2 vertices with degree 1 
and n — 2 vertices with degree 2. For highly connected clusters, such as cubic clusters 
in d > 1 dimensions (see Remark 2.5 below for their definition), the bound (2.6) is 
very crude as compared to (2.5). Though, (2.6) will suffice for our purpose. 
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The next lemma provides a Faber-Krahn inequality on graphs. In contrast to 
Cheeger inequalities, such estimates for graphs have not been known for a long time, 
see [13] for a detailed exposition. Lemma 2.4 adapts a result from [13], which is 
proven there for more general graphs, to the type of graph Laplacians we use here. 

Lemma 2.4 (Faber-Krahn inequality). Let := (Y, <S) be a connected finite sub- 
graph ofh d with \y\~^2 vertices. Then its lowest Pseudo-Dirichlet eigenvalue obeys 

where /ifk G]0, oo[ is a constant that depends only on the spatial dimension d. 

Remark 2.5. The Faber-Krahn inequality produces asymptotically the correct In- 
dependence as \Y\ — > oo, if, for example, ^ is a cubic cluster J3i, that is, £li = 
for some finite cube Ai C Z d with \ A\\ = l d points, i.e. edges of length i - 1 6 N. 

Proof (of Lemma 2.4). We reduce the assertion to a particular case of Prop. 7.1 and 
Cor. 6.4 in [13] by choosing the weighted graph in Prop. 7.1 as L d with unit weights 
on all bonds - note that these results in [13] extend to d = 1. Given any A C Z d , this 
yields the inequality 

M4:= ^ ^^Z® (2-8) 
V ; o^ v ec„(A) 2d{y,<p) 2\A\ 2 / d 

The constant (3d G]0, oo[ is the isoperimetric constant of Cor. 6.4 in [13], which is 
independent of A C Z d . Moreover, co(A) stands for the £ 2 (Z d )-subspace of real- 
valued sequences with support in A. In order to check that the above definition of 
Ai(yl) matches the one in [13], we refer to Sect. 5.5 of that paper. The claim of the 
lemma follows from the estimate 

Eg ) (9)>2d\ 1 (r), (2.9) 

which we prove now. To this end we observe that the adjacency operator has only non- 
negative matrix elements (5 X , A(%?)5 y ) for all i,]/ 6 f and that A(^) is irreducible 
on £ 2 (y) due to the connectedness of . Hence, it follows from the Perron-Frobenius 
theorem, see e.g. Thm. 2.1.4 in [3], that the eigenvector of A^S) = 2di - A{<£) 

corresponding to the non-degenerate smallest eigenvalue ('£) can be chosen to 
have non-negative entries. This implies 

inf ^f 5 ^- (2.10) 

0(a;)>O for all xe'f 

The proof of (2.9) is completed by noting that 

(<t>,A^)<t>) = ]T [^)-0(y)] 2 +^[2rf-^(x)](0(x)) 2 

[x,y\eS x^'f 
[x,y]eE d [x,y}e¥. d \g 

^(^Z\ N (L d )^) (2.11) 
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for all 4> as in (2. 10), where <p(x) := { ^ ' ^ J \ r □ 

As a last ingredient for the proof of Theorem 1 . 14 we need some simple upper es- 
timates on e£ ifS) for special types of clusters. These estimates are obtained from the 
minmax-principle. For our purpose it is only important that they reproduce the asymp- 
totically correct functional dependence on the number of vertices for large clusters. 

Lemma 2.6. (i) The lowest non-zero Neumann eigenvalue for a linear cluster 

Jz?„, which was defined in Remark 2.3(ii), obeys 

4 1} (^nK^ (2.12) 

for allneKn^ 2. 

(ii) The lowest Dirichlet eigenvalue for a cubic cluster i?;, which was defined 
in Remark 2.5, obeys 

4 1} (^K^ (2-13) 

for all I G N, I ^ 2. 

Proof. (i) The minmax-principle yields the upper estimate 

n—lr \2 



V-V71— 1 / 



ii i 



4 1} (J2f n ) ^ 2 - (2. 14) 

for every (m, . . . , u n ) G W l subject to the orthogonality constraint Y^j=i u i = 0- 
Choosing Uj := — j + (n + l)/2 for j G {1, . . . , n}, which is in accordance with the 
orthogonality constraint, proves part (i). 

(ii) Appealing to the minmax-principle with a trial "function" that factorises 
with respect to the d Cartesian directions, gives 

— Jr j , — — (2.i5) 

for all (tti, . . . S M'. Now, we choose w 3 := -|j - (I + 1)/2| + (Z - l)/2 for all 
j G {1, ...,/} so that ui = u; = 0. If I ^ 3 is odd, an explicit calculation shows 

while for an even integer Z ^ 2 it yields 



and the lemma is proven. □ 
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The next two lemmas provide the key estimates for Theorem 1.14. While the lower 
bounds in Lemma 2.9 hold for all p e]0, 1[, the upper bounds in Lemma 2.7 are 
restricted to the non-percolating phase. 

Lemma 2.7 (Upper bounds). Let d e N and consider p e]0,p c [. Then there exist 
constants ot^ , af^ G]0, oo[ such that 

N N (E) - N N (0) < exp{-a+£;- 1/2 } for all E e]0,4rf] , (2.18) 
Nu(E) N b (E) < exp{-a+E- d/2 } for all E e]0,2d[. (2.19) 

Remark 2.8. It is only the right inequality in (2.19) whose validity is restricted to E 6 
]0,2d[. The proof below will show that N D (E) ^ cxp{-a+ir d / 2 } holds for all 

Ee]0,4d[. 

Lemma 2.9 (Lower bounds). Let d e N and consider p g]0, 1[. Then there exist 
constants G]0, oo[ such that for every E e]0, 4<i] one was 

N N (E) - N N (0) > exp{-aN^" 1/2 } , (2.20) 
iV B (£) ^ iV D (S) > exp{-a D S- d / 2 } . (2.21) 

froo/ (of Theorem 1.14). Due to the symmetries (1.21), it suffices to prove the as- 
serted Lifshits tails at the lower spectral edge. These follow from the estimates in 
Lemma 2.7 and Lemma 2.9, because after taking appropriate logarithms, the respec- 
tive bounds coincide in the limit E J. 0. □ 

Proof (of Lemma 2.1). Fix E e]0,4d] subject to E < 7x if X e {D,D}. The con- 
stants 7x were defined in Remark 1.5(iv). Definition 1.9 of the integrated density of 
states implies 

N X (E) - N x (0) = f P(dco) (J , [0(E - A^) - P^] Jo) , (2-22) 
Jn 

where Px ■= 0(—Ax) denotes the (random) projector onto the null space of Ax- Due 
to our assumptions on E and the block-diagonal form ( 1 . 14) of Z\ x , the right-hand side 
of (2.22) is only different from zero if the origin is part of a cluster with at least two 
vertices. Let us call this event f2o and the corresponding cluster :— (Yq^ 1 , 4"' ) 
for all cj E J?o- Hence, we obtain 

N X (E) - JVx(O) = / F(dw) (Jo, [0(E - A x {^ u) )) - Px(^o M )]<5o) 
Jn 

^ f P(dcj) 0(E - 4 1) (^ ) ))( J , [1 - Px(^o M )]^o) 
J n 

a{w€flb:£>4V)}' (2.23) 

Before we make a distinction of the two cases X = N and X £ {D, D} in order to 
apply the Cheeger, respectively the Faber-Krahn inequality, we recall from Proposi- 
tion 1.2 (and Remark 1.3(h) for the case d — 1) that the cluster is finite for 
P-almost all w G Qq, since we assume p < p c . 
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Neumann case. Applying the weakened version (2.6) of Cheeger's inequality to 
(2.23), yields the claim 

N N (E)-N N (0) ^ P{co e n : |r (w) | ^ 1/idE) 1 ' 2 } < exp{-<r ^(p)^-^} . 

(2.24) 

The second inequality in (2.24) reflects the exponential decay of the cluster-size dis- 
tribution in the non-percolating phase, see Thm. 6.75 in [16]. Here, ((p) > is some 
finite constant for every p G]0, p c [, which depends only on d. Formally, Thm. 6.75 in 
[16] does not cover the one-dimensional situation d = 1 . But for d = 1 the exponential 
decay of the cluster-size distribution follows from elementary combinatorics. 

(Pseudo-) Dirichlet case. The inequalities (1.18), (2.23), the Faber-Krahn in- 
equality of Lemma 2.4 and the exponential decay of the cluster-size distribution yield 

N D (E) < N b (E) < F{uj e f2o : l^l > (h FK /E) d ^} < C xp{-h d F £c(p)E- d/2 } . 

(2.25) 

Here we have used that 7Y D (0) = (0) = 0, see Remark 1 . 15(ii). □ 

Proof (of Lemma 2.9). Lemma 1.12, the isotony of the right-hand side of (1.22) in E, 
the right-continuity of Nx and (1.25) imply that 

N X (E) - ^Vx(O) ^ limsup ±- trace, 2(/1) [0(E - A%\) - p^\] (2.26) 

A\l d \ A \ 

for all E > and all u> € fl 1 . Concerning the limit in (2.26), we think of a sequence of 
expanding cubes that are centred at the origin, cf. Remark 1.1 3(i). Due to the block- 
diagonal form (1.14) of Z\ x , which continues to hold for with respect to the 

decomposition of ^ into clusters V^J := j £ {!,•••, J^}, we 

get 

N X (E)-N X (0) > limsup -j- V trace, 2 M [G{E- A x {^)) . 

A\TA \ A \ ^ y A 'i 1 ,J ,J 

(2.27) 

At this point we make again a distinction of the cases X = N and X e {D, D}. 

Neumann case. Let £„ denote the set of all linear clusters with n ^ 2 vertices 
in L d and let Xc n be the characteristic function of this set of graphs. A crude lower 
bound on the j-sum in (2.27) for X = N results from discarding all branched clusters, 
i.e. those which are not linear, 

•'a oo 

££Xc B (*t?) t^e p{<)) [0{E-A N {^])) -P N (<J)] 

j=l n=2 

oo -'a 

^££^(^)0(^-4 1) (^S)) 

n=2 j=l 
oo 2 

£ £ 0(£ - 12/n 2 ) ^ X £ „ (tfg) . (2.28) 

n=2 j=l 
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The first inequality in (2.28) follows from restricting the trace to the spectral subspace 
corresponding to E^p { c to^ ), the second inequality follows from the variational upper 
bound in Lemma 2.6(i). Thus, (2.27) and (2.28) yield 

oo 

N N (E) - JV N (0) > limsup V0(£- 12/n 2 ) n" 1 L^(yl) (2.29) 

with 

L^(A) : = |S]E X ^{<]) = ^ A : <4"\x) € (2.30) 

being the number density of points in A that are vertices of a cluster of type £„. Here, 
(x) denotes the cluster of that contains x e A. For a given n e N, n ^ 2, 
and a sufficiently large bounded cube 7l C Z d with {A] 1 / 4 ^ 2n + 1, let us also define 
the number density 

L { n\A) :=^-\{x&A: min min \x v - y v \ > n+ 1 and ^(z) £ £„)| 

(2.31) 

of vertices which are, in addition, sufficiently far away from the boundary of A. 
Clearly, one has 

\\m[L^{A)-L^{A)]=Q (2.32) 

for all u <E fl and all n E N, n ^ 2, since the difference in the two quantities results 
from a surface effect. The vertices that count for (A) are so far away from the 
boundary of A that the clusters they belong to cannot grow when enlarging A. Hence, 

\A 1 UA 2 \L^(A 1 UA 2 ) > \A 1 \L^(A 1 ) + \A 2 \L^(A 2 ) (2.33) 

holds for all Ai,A 2 C 7L d provided A\V\A 2 = 0. Thus, L n (A) defines a superergodic 
process and the Ackoglu-Krengel superergodic theorem, see e.g. Thm. VI. 1.7 in [10], 
and (2.32) imply 

lim L^(A) = sup / P(dw') L^'\A) = ¥{lu' 6 Q : %T (a/) e £„} (2.34) 

for all n G N, n ^ 2, and P-almost all co € O. The event f2 and the random cluster 
were defined above Eq. (2.23). 

Now, we neglect all terms in the n-sum in (2.29) except for the one which corre- 
sponds to the biggest integer n(E) obeying n(E) < (12/E) 1 ' 2 + 1. From this we 
conclude together with (2.34) that 

N N (E) - 7V N (0) ^ [n(E)}- 1 P{uj e Q : ^ G &n(E)} ■ (2.35) 

Elementary combinatorics shows that the probability on the right-hand side of (2.35) is 
bounded below by cxp{— n(E) f(p)}, where f(p) G]0, oo[ is a constant that depends 
only on d for a given p e]0, 1[. This leads to the estimate 

N N (E)-N N (0)^ (12/ e ^ + 1 expj-12 1 ^ f{p)E~^}, (2.36) 
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which can be cast into the form (2.20) for E e]0, Ad]. 

(Pseudo-) Dirichlet case. This case parallels exactly the previous one, except that 
here we retain cubic clusters instead of linear clusters. Let £2/ denote the set of all 
cubic clusters in L d with l d vertices, i.e. edges of length I — 1 e N, let Xq, be the 
characteristic function of this set of graphs and define the number density (A) := 

(l d /\A\) X)j=i -^Qi C^A^j) °f points in A that are vertices of such a cubic cluster (when 
restricted to A). Now, the role of £„ in the previous case will be played by 0;. Hence, 
the analogue of (2.29) reads 

oo 

N b (E) ^N U (E) ^limsu P V<9(£-27d/Z 2 )r d Q ; M (yl), (2.37) 

where we have used Lemma 2.6(ii) instead of Lemma 2.6(i). The very same arguments 
that led to (2.34) imply in the present context 

lim Q, M (A) = ¥{lo' G A) : c $t' ) e &i } (2.38) 

for all I G N, I ^ 2, and P-almost all u G fl. By neglecting all terms in the ^-sum 
in (2.37) except for the one which corresponds to the biggest integer 1(E) obeying 
1(E) < (27d/E) x / 2 + 1, we conclude with (2.38) that 

N b (E) > N U (E) > [l(E)]- d P{uj e f2 : ^ e Qi(^)} • (2-39) 

Again, there is an elementary combinatorial lower bound exp{— [l(E)] d g(p)} for the 
probability in (2.39), where g(p) g]0, oo[ is a constant that depends only on d for a 
given p g]0, 1[. So we arrive at 

7V B (i?) ^ 7Y D (S) ^ [(27d/£;) 1 / 2 + l]- d cxp{-[(27d/£;) 1 / 2 + l] d g( P )} , (2.40) 

which can be cast into the form (2.21) for E g]0, Ad] . □ 
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